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Abstract 

We develop a covariant approach to describe the low-lying scalar, pseu- 
doscalar, vector and axialvector mesons as quark-antiquark bound states. 
This approach is based on an effective interaction modeling of the non- 
perturbative structure of the gluon propagator that enters the quark 
Schwinger-Dyson and meson Bethe-Salpeter equations. We consistently treat 
these integral equations by precisely implementing the quark propagator func- 
tions that solve the Schwinger-Dyson equations into the Bethe-Salpeter equa- 
tions in the relevant kinematical region. We extract the meson masses and 
compute the pion and kaon decay constants. We obtain a quantitatively 
correct description for pions, kaons and vector mesons while the calculated 
spectra of scalar and axialvector mesons suggest that their structure is more 
complex than being quark-antiquark bound states. 
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I. INTRODUCTION 



Recently, the scalar mesons have attracted a lot of interest as the reanalysis of the pseu- 
doscalar meson scattering data indicated the existence of a flavor SU(3) nonet in this chan- 
nel [Q . It is therefore desirable to gain deeper understanding of the constituent structure of 
the scalar mesons together with a comprehensive description of the meson states in the other 
spin-parity channels. The ultimate goal would be to understand all low-lying meson states 
and resonances as non-perturbative bound states in Quantum Chromo Dynamics (QCD). 

A relativistic framework for analyzing mesons as composite objects is provided by the 
Bethe-Salpeter equations that extract poles in the quark-antiquark scattering kernel 0. 
The attraction needed to bind quarks and antiquarks emerges from dressed multiple gluon 
exchange. Thus the essential ingredients to these equations are the quark and gluon prop- 
agators as well as the quark-gluon vertex. In addition, these ra-point Green's functions are 
related by their Schwinger-Dyson equations which are part of an infinite tower of non-linear 
integral equations. There has been some progress in the understanding of the infrared be- 
havior of the gluon propagator from recent Yang-Mills lattice calculations as well as from 
studies of the coupled system of gluon and ghost Schwinger-Dyson equations Never- 
theless, for phenomenological applications the frequently adopted strategy is to model the 
gluon propagator as well as the quark-gluon vertex and consistently derive the quark propa- 
gator from its Schwinger-Dyson equation because it facilitates the continuation to complex 
Euclidean momenta. 

These types of calculations have a long history, for reviews see Refs. 0,0. Early versions 
adopted pointlike gluon propagators in coordinate space that eventually lead to Nambu- 
Jona-Lasinio (NJL) type models pointlike propagators in momentum space were also 
considered [§. These models are particularly simple because either solving the Schwinger- 
Dyson equation yields a free quark propagator or the Bethe-Salpeter integral equations 
reduce to algebraic equations. The main target particularly of the NJL-model studies have 
been the pseudoscalar mesons. It turned out that they can be adequately described once 
the important feature of dynamical chiral symmetry breaking is incorporated, i.e. the in- 
teraction is strong enough so that the resulting quark propagator signals a non-vanishing 
quark condensate reflected by a non-zero constituent quark mass. Then the pseudoscalar 
mesons can be understood as the would-be Goldstone bosons of chiral symmetry breaking. 
However, as the Bethe-Salpeter equations involve the dressed (constituent) quark propaga- 
tors, binding can only be achieved kinematically and meson states with masses larger than 
twice the constituent quark mass cannot be described consistently. For that reason, model 
gluon propagators have been developed that yield quark propagators without poles for real 
momenta as an attempt to include the confinement phenomena PJT(][]. Again these studies 
focused on pseudoscalar mesons 0)111] while a comprehensive investigation for the scalar, 
pseudoscalar, vector and axialvector mesons has not been carried out so far. Other stud- 
ies P,p!0| made contact with perturbative QCD by considering a model gluon propagator that 
matches the pertinent anomalous dimension. This contribution has negligible effect on the 
meson properties, but its inclusion makes cumbersome the extraction of the solutions to the 
Schwinger-Dyson equations for the large time-like momenta that enter the Bethe-Salpeter 
equations with meson states other than the pseudoscalars. We consider this attempt in- 
teresting, but as we do not want to focus on the ultraviolet properties of mesons here, we 
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regard it an unnecessary technical complication. Therefore we focus rather on establishing 
a model as simple as possible that we consider a pertinent starting point to study the struc- 
ture and properties of low-lying mesons in a relativistic framework. Our model interaction is 
parameterized in form of a non-trivial gluon propagator that contains sufficient strength to 
cause dynamical chiral symmetry breaking. For technical reasons it turns out that a Gaus- 
sian shape function for the propagator in momentum space is most suitable. Essentially we 
consider this model propagator as an effective interaction that relativistically describes the 
binding of quarks and antiquarks to mesons. Furthermore, we take the quark-gluon vertex 
function to be the tree level one since this procedure provides a framework that is consistent 
with chiral symmetry when the ladder approximation for the Bethe-Salpeter equation is 
employed For approaches going beyond ladder approximation see e.g. Ref. [|T^] . 

This paper is organized as follows. In Section II we will introduce the effective interaction 
and solve the Schwinger-Dyson equation for the quarks. We will put particular emphasis on 
the analytic continuation of the resulting quark propagator to time-like momenta that enter 
the Bethe-Salpeter equations. We will discuss the structure of the Bethe-Salpeter equations 
and present the solutions in Section III. In Section IV we will present the numerical results 
in the sector of the three light quarks (up, down and strange) and compare them to the 
empirical data. In Section V we will conclude and suggest a possible extension of the 
current approach in particular with regard to the possibility that the scalar meson might 
have to be considered as two-quark - two-antiquark bound states We devote an 

appendix to discuss the numerical stability of our results. 



II. THE QUARK SCHWINGER-DYSON EQUATION 

The purpose of this section is twofold. First we will explain our model gluon propagator 
whose dressing parameterizes the effective interaction and then solve the corresponding 
Schwinger-Dyson equation for the quark propagator. 

As already discussed in the introduction we take a Gaussian form for dressing the model 
gluon propagator. This follows the work of Ref. [|l^ with the exception that we omit 
the logarithmic tail that matches perturbative QCD because we are not interested in the 
ultraviolet properties of mesons. We therefore write. 



where ix,v are Lorentz indices, t^j^{q) is the transverse momentum projector and a, 6 label 
color. While the prefactors in eq. are chosen to make subsequent equations more con- 
cise, D and uj are dimensionful parameters that we will determine from fitting empirical 
data. The coefficient D sets the strength of the interaction and uj is the value at which the 
scalar function in the parameterization is maximal. Hence uj sets the interaction scale. The 
dressed gluon propagator is supposed to represent a sensible hadron model and hence 
one can envisage that uj will have a value of several hundred MeV. 

We interpret the effective interaction (0) as the propagator (in Landau gauge) of a gluon 
that gets absorbed and emitted by the quarks that eventually get bound to form mesons. 
To completely define the interaction, we need to parameterize the quark-gluon coupling. 
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To establish chiral symmetry we apply the rainbow-ladder approximation to the system of 
Schwinger-Dyson and Bethe-Salpeter equations. This implies that the quark-gluon coupling 
is given by the tree level interaction vertex, ig^y^j^^, where A" is a Gell-Mann matrix acting 
in color space. Note, that we have already included the coupling constant g in the definition 
of the effective interaction (|T]). 

Then the Schwinger-Dyson equation for the (inverse) quark propagator becomes [d^k = 



S~-\p) = ^^ + mo + J d'ky,S{kh.g'^^Gf,{k-p) 



(2) 



where mg is the current mass of the considered quark. This contribution represents the only 
explicit distinction between quarks of different flavors. Of course, its effects will implicitly 
propagate through the whole calculation. However, for notational simplicity we will continue 
to suppress flavor labels. A suitable parameterization of the quark propagator is inspired 
by the form of a free fermion propagator 
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(3) 



In solving the Schwinger-Dyson equation (H) we have to flnd the scalar functions A{p'^) and 
B{p'^). It is also very instructive to deflne a mass function via M{p'^) = B{p'^) / A{p^). In 
particular M(p^ = 0) plays the role of a constituent quark mass and a large value thereof 
signals dynamical chiral symmetry breaking. 

As usual we work in Euclidean space with Hermitian Dirac matrices^ that obey {7^, 71,} = 
25^j, and 75 = —71727374. Inserting the effective interaction (|l]) and performing the standard 
trace algebra, we then deduce the following coupled equations for the propagator functions 



A{x) 



D 

1 + ^ 



dyyA{y) 



uj^Jo (yA^y) + B^{y)) 



X— / dzy/l — 



3^ 



1 + ^ 

X 



xyz~ -yz 



exp 



X + y — 2^Jxy z 



(4) 



B{x) 



D 

mo + — 



dyyB{y) 



ij^Jo (yA^iy) + B^iy)) 



: — J dz\/l — z'^ [x + y — 2^xy z\ exp | 



X ^ y — 2Jxyz 



(5) 



where the four dimensional integral measure has been expanded such that x = p^^y = k"^ 
and z = p-kj y/jPW. 

In a flrst step we solve eqs. @ and (|) for spacehke momenta, i.e. for real positive x. 
Then we observe that the integrals on the RHS of these equations only involve the prop- 
agator functions at real arguments y and we can use them to numerically compute the 



^This can be related to the standard Minkowski space Dirac matrices via 74 



-7f. 



4 



propagator functions for arbitrary complex x. At first sight, it appears that A{x) and B{x) 
could not be consistently continued because the cut along the negative x-axis (associated 
with ^/x) would yield different results when continuing in the upper or the lower half- 
plane and it would be impossible to resolve the ambiguity in ^/x when continuing 
X Fortunately, this is not an obstacle. By expanding the exponential functions in 
eqs. (§) and it becomes obvious that all terms that contain ^/x are odd in z and thus 
vanish when integrating over this angular variable. Thus we are free to choose either of the 
two signs above. For definiteness we work in the upper half-plane with ^/x —>■ along 
the negative half-line. 

The observation that it is sufficient to know the propagator functions along the real 
positive axis to compute them for all complex arguments of interest is a major point of this 
present study. It is an especially important issue to the study of the Bethe-Salpeter equation, 
since the quark propagator must be sampled in the complex plane of Euclidean momenta. In 
other words, the quark functions are calculated at precisely the momenta for which they are 
used, with no fitting functions, interpolation or extrapolation. However, in order to do so 
numerically, the quark propagator functions must be known very accurately along the real 
positive axis. This is achieved by noting that in eqs. (^) and (^) the angular {z) integrals 
can be done analytically for real and positive x. This property of the angular integrals is a 
particular feature of the Gaussian dressing function (0). The resulting equations read 



{yA\y) + B\y)) "I uj^ 




dyyB{y) r x + y 



u( \ ^ n r dyyB{y) 



,2 




{yA^{y) + B^{y)) "I 

(6) 

where J„ are modified Bessel functions. These equations are one-dimensional coupled non- 
linear integral equations which can be straightforwardly solved numerically with a high 
degree of precision. The solutions for real (Euclidean) momenta are subsequently substituted 
into eqs. (H) and (B) to yield the propagator functions for complex momenta. 



III. THE BETHE-SALPETER EQUATION 



Having obtained the quark propagators in the complex plane from the Schwinger-Dyson 
equations we have collected all ingredients for the Bethe-Salpeter integral equations. They 
will ultimately yield the quark meson vertex functions, F, that describe mesons as bound 
quark-ant iquark pairs. Strictly speaking, the Bethe-Salpeter equation is an eigenvalue prob- 
lem, valid only at the resonance pole = — M^, where M is the mass of the resonance. It 
is derived from considering the four-point quark Green's function that involves exchanges 
of resonance mesons. Such an exchange is characterized by a pole in that four-point func- 
tion, and the homogeneous Bethe-Salpeter equation below (^ determines the position of 
that pole. All other regular terms in the vicinity of this pole are neglected. Demanding 
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furthermore that the residue of this pole is unity yields the normalization condition ([T^ ) for 
the vertex functions. 



A. Bethe Salpeter vertex functions 



The vertex functions resulting from the Bethe-Salpeter equation are characterized by 
three momenta out of which only two are linearly independent due to momentum conser- 
vation at the vertex. If we denote the meson momentum P and the momentum of the 
incoming quark p + C,P then the momentum of the outgoing quark (= incoming antiquark) 
is p+ — 1)P. This suggests to label the vertex functions by p and P: T{p, P) . We have also 
introduced the arbitrary momentum partition parameter ^ e [0, 1]. Due to strict relativistic 
covariance the results for physical observables do not depend on ^. Unfortunately we will 
have to assume approximations within the full numerical computation (see Section 4) that 
violate covariance to some extent. We will study the ^-dependence of our results and verify 
that within a wide range ^-independence is maintained, see Appendix A. This will represent 
an a posteriori validation for the relativistic covariance of our computations. 

We now turn to the main target of our studies, the Bethe-Salpeter integral equations for 
the vertex function r{p,P) in ladder approximation P,p|,p!5|: 

T{p; P) = -y d'k [7. S{k + ^P) T{k; P) S{k + - 1)P) 7^] g'G.^k - p) . (7) 

Here we have factorized the color factors in the effective interaction, G'^^^{q) = 5"'^G^p{q) and 
performed the corresponding trace. The flavor content of the meson is not made explicit in 
eq. as we have suppressed the flavor labels in the quark propagators. It is understood 
that the two propagators in eq. (|^) are taken such as to account for the flavor quantum 
numbers of the considered meson. In the model that we will consider, the up and down 
quarks will be assumed to have equal current masses (mo in eq. (^) and thus also identical 
propagator functions A{x) and B{x). For the light quarks, which should give rise to the 
familiar SU{3)f nonet, we are thus left with three representatives of each of the multiplets 
that are distinguished by their isospin number, / = 0, |, 1 . We must also specify the meson 
angular momentum and parity. This is reflected by the Dirac and Lorentz decomposition 
of the meson vertex functions. This decomposition is known in the literature and here we 
follow Ref. [0. For the pseudoscalar channel (J^ = 0~) we take 



r(^)(p; P) = 75 rlf\p; P) - irvrip; p) - ^i>T^{p■, P) 



The decomposition for a scalar (J^ = 0^) meson reads 

r(^)(p; P) = T^^\p; P) - ^rT['\p■, P) - ^i^Pf (p; P) - [fj] if ^(p; P). 



(8) 



(9) 



The vector {J^ = 



channel involves eight scalar functions 



ri^) (p; p) + rTY> {p- p) - i^r^^ ' {p- p) + 1 [f , i>] T'^ > {p- p) 



+ 
+ 



p2 

P^^P-P 
p2 

P^P-P 
p2 



V^P{p-P) + 2i]f>T'i\p-P) 



(V), 



(10) 



rf ) {p-P)+ i]f>T\' \p-P)- ii>T),' \p;P) + [f , ^] >{p-P) 



iv). 
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In the axialvector channel we have two modes that are distinguished by their charge conju- 

For = : 



gation properties []16 



7m 



zTi''\p; P) + rvrip; P) - ^T^'ip; P) + i [f,i>] Vf\p- P) 



the decomposition is 



+75 



Pa* 



p2 
p2 



(A), 



v[^\p-P) + 2i]/^Vf\p-P) 



{A), 



(11) 



while the J^*" 



mode is decomposed as 



^i^\p:P) 



75 



p2 



t[^\p- P) + ifT'C\p- P) - tiyVC'ip- P) + [f.i)] rr(p; P) 



{A), 



(12) 



In what follows we will omit the superscripts that label the spin and parity channels because 
these channels do not mix and there should hence be no confusion. 

The primary object of this paper is to extract the bound state masses for the various 
flavor combinations and angular momentum channels.0 The corresponding projection results 
in sets of coupled equations for the Fj. After carrying out two of the three angular integrals 
analytically we are left with functions of the squared momenta p^ and P^ as well as the 
angle between p and P: z = p-P/ \/p^P'^. The 2;-dependence is analyzed by an expansion in 
Chebyshev polynomials 



Up;P) 



E 

k=0 



\fc-pfc 



(13) 



Since the form an orthonormal set, we can project the equations for the Dirac components 
onto r^(p^; P^). Finally, the fc^-integral in the Bethe-Salpeter equation (^ is implemented 
numerically as a matrix equation for the unknown r^(p^;P^), p'j being the discrete values 
of the momentum squared. The kernel, K of that matrix parametrically depends on the 
meson momentum P^. We solve that matrix equation as an eigenvalue problem by tuning 
the meson momentum to P^ = — M^, such that Det(l — K) = 0. This yields the desired 
meson mass M. 

For practical reasons we need to truncate the expansion (|13D at a certain order, Nch- 
As noted before this violates relativistic covariance. Fortunately the momentum partition 
dependence provides a measure for the degree at which covariance is violated. We can 
mitigate this degree by increasing Net- Numerically we furthermore scan the momentum 
partition dependence of the computed meson masses to a posteriori verify covariance for our 
solutions to the Bethe-Salpeter equations, see Appendix A. Unless otherwise stated we use 
Nch = 4. 



^ The traceology that is involved to project the Bethe-Salpeter equation (^) onto the Dirac com- 
ponents Fj is lengthy, but straightforwardly performed using the algebraic manipulation package 
FORM ^ 



7 



B. Pseudoscalar leptonic decay constant 



The solution to the Bethe-Salpeter equation not only yields the meson masses but also 
the meson quark vertex functions that can be used to compute meson properties. Here 
we will focus on the pseudoscalar decay constants and fx- In order to calculate these, 
we first have to normalize the vertex functions r(p;P). The Bethe-Salpeter equation is a 
homogeneous equation, and thus needs an additional normalization condition. As mentioned 
previously, that condition is obtained from demanding the pole in the four-quark Green's 
function to be unity. For equal momentum partitioning, {i.e. for ^ = 1/2 only) it reads fl^ 

2P^ = 3| d'kTil^Tik, -P)^^^t^^Tik,P)Sik- P/2) 

+T{k, -P)S{k + P/2)r(fc, p f^%~p^^^^ ] (14) 

where the trace is over Dirac matrices. The conjugate vertex function F is defined as 

T{p,-P)=CT^{-p,-P)C-\ (15) 

where use is made of the charge conjugation matrix, C = —7274. The quark propagator 
derivatives are calculated by differentiating the quark Schwinger-Dyson equations (H) and (^ 
analytically and then numerically integrating the corresponding expressions. 

The decay constants are finally obtained from the coupling of the axial current to the 
quark loop |jl5| 



f = jpj d'kTT{r{k,-P)Sik + P/2)^,rS{k-P/2)}, (16) 

where again the trace is over the Dirac matrices and the vertex functions are normalized 
according to eq. (p!^. Again we have not made explicit the flavor quantum numbers. 



IV. NUMERICAL RESULTS AND DISCUSSION 

We are now in the position to present our numerical results to the system of Schwinger- 
Dyson and Bethe-Salpeter equations. Even though the meson results are the primary target 
of our investigation it is fruitful to first consider the quark propagator functions obtained 
from the Schwinger-Dyson equations, as they enter the kernel of the Bethe-Salpeter equa- 
tions. 



A. Quark propagator functions 

In Fig. |1] we show the quark propagator functions v4(p^), -B(p^) and M(p^) = B{p'^)/A{p'^) 
along the positive, real spacelike axis, p"^ > 0. This is the numerical solution to the coupled 
equations (|^), which we emphasize is the basis for the quark solutions for complex momenta. 
The solution clearly shows that dynamical chiral symmetry breaking is occurring: the mass 
function M{p^) attains a sizable non-zero value, even in the case that the bare quark mass 
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FIG. 1. Plot of the spacelike chiral quark functions (as a function of the momentum squared). 
The parameters are uj = O.SGeV, D = 16.0GeV~^. Left panel: mo = 0, right panel: mo = 0.115 
GeV. All units are GeV. 



UJ 


D 


mo 


M{p^ = 0) 


0.4 


45.0 





0.520 






5 X 10^3 


0.531 






0.12 


0.709 


0.45 


25.0 





0.450 






5 X 10^3 


0.462 






0.12 


0.657 


0.5 


16.0 





0.397 






5 X IQ-^ 


0.413 






0.115 


0.622 



TABLE L The variation of the quark mass function evaluated at zero momentum squared. 
All units are GeV. 

mo is zero. This phenomenon has been extensively studied (see for example Refs. 
and is an example of genuinely non-perturbative behavior as dynamical mass generation 
cannot occur at any order in perturbation theory. Recall that the effective interaction (|l]) 
that enters the Schwinger-Dyson equations does not contain the perturbative UV behavior, 
rather it has an exponential damping at high momenta. This is manifested in the quark 
propagator functions as a sharp transition from the low momentum behavior to the bare 
values in the high momentum region. This transition occurs at about IGeV. 

The quark mass function M (p^) is also a useful quantity to discuss the changes of the 
quark propagator with the input parameters. Table | shows M(p^ = 0) for those parameter 
sets oi uo, D and itlq that we will later use to compute the meson masses. As already 
mentioned the quark mass function indicates the extent to which dynamical chiral symmetry 
breaking occurs and plays the role of the mass of the constituent quarks within hadrons. 
We see that the quark mass function decreases as the effective scale uj increases and the 
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FIG. 2. Plot of the timelike quark mass function squared M'^{p^). The parameters are 
uj = 0.5, D = 16.0. Left panel: mo = 0, right panel: mo = 0.115. All units are GeV. 



strength D decreases. 

Since the main thrust of the paper is the use of the quark propagator functions eval- 
uated at precisely the complex momenta squared for which they enter the Bethe-Salpeter 
equation, it is important to study the behavior of the mass function at least along the time- 
like axis, i.e. for negative real momenta squared. In Fig. ^ we present typical results for 
the mass function squared M^(p^) in that regime. Although the effective interaction (|I]) is 
simple, the timelike functions are distinctly non-trivial. Already from Fig. |^ we observe a 
significant variation with the current quark mass. This indicates that the specific forms of 
timelike quark propagator functions possess pronounced model and parameter dependences. 
Investigating the model dependence requires changing the effective interaction (|I]) and is not 
subject of the present study. For the parameter dependence, we fortunately find (see be- 
low) that it does not significantly effect the model predictions for the meson masses. There 
are divergences whenever A and B change sign at different p^. Note that there are several 
particle-like poles (p^ = — M^(p^)) where the quark goes on-shell. This is an important 
observation, since it is not at all clear whether a confined quark can be on-shell. Indeed, 
one might be tempted to conclude that with on-shell singularities and a subsequent lack of 
confinement, the model will not be able to give meson masses at all. As we will see, this is 
fortunately not the case. 



B. Mesons 



Turning to the mesons, we note that there are four model parameters, uj,D, niu and rus 
that we first have to fit to empirical data. To this end, we initially choose the pseudoscalar 
meson observables M^, Mk and /,r- Then one parameter remains unconstrained by the pseu- 
doscalar sector alone. However, the condition that the quark propagator function reflects 
dynamical chiral symmetry breaking, i.e. that M(p^ = 0) ~ O.SGeV leaves only a small 
window for the remaining choice, c/. table |. All other masses and decay constants are sub- 
sequently model predictions. The resulting model parameter and the predicted kaon decay 
constant fx, that is unexpectedly well reproduced, are shown in Table ||. The subsequently 
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D 


ruu 


ms 






Mk 


fK 


0.40 


45.0 


5 X 10-3 


0.120 


0.135 


0.131 


0.496 


0.164 


0.45 


25.0 


5 X 10-3 


0.120 


0.135 


0.131 


0.496 


0.163 


0.50 


16.0 


5 X 10-3 


0.115 


0.137 


0.133 


0.492 


0.164 


experiment [1^] 


0.135 


0.131 


0.498 


0.160 



TABLE II. Parameter sets used and fit results for the pseudoscalar mesons. M^^, /^r and Mk 
are used as input, fx is predicted. All units are GeV. 



UJ 


D 


niu 


nis 


Mp 


Mk- 




0.40 


45.0 


5 X 10-3 


0.120 


0.748 


0.939 


1.072 


0.45 


25.0 


5 X 10-3 


0.120 


0.746 


0.936 


1.070 


0.50 


16.0 


5 X 10-3 


0.115 


0.758 


0.946 


1.078 


experiment [18 


0.770 


0.892 


1.020 



TABLE III. Results for the vector mesons. All units are GeV. 



predicted meson masses are shown in Tables |ITI| - |V]] . In all cases we have the inequalities 



Muu < Mus < Mssi where the subscript labels the flavor content. These relations just reflect 
the quark-ant iquark picture that is implicit in the present Bethe-Salpeter approach. 

Obviously both the pseudoscalar (table |l|) and vector mesons (table pT|) can be very 
well described within our model with the choice uj ~ 0.5GeV. Our results agree with a 
previous analysis of the vector mesons based on an effective interaction which included the 
perturbative type term ||TD[. This shows that such terms do not have a large effect on the 
meson masses, at least for the pseudoscalar and vector cases. Indeed, in the context of 
low-energy meson phenomenology we conclude that the logarithmic tail, and its associated 
renormalization represent an unnecessary obfuscation. 

The situation for the scalar mesons (table [TVD is not quite that clear. To begin with the 
particle data group [0 does not provide a clear picture in this channel but only quotes a 
wide range for the mass of the lowest scalar (0.4 - 1.2 GeV). More detailed studies of the 
pseudoscalar scattering amplitudes revealed that the assignment of the scalar meson nonet 
is not at all established |jl| . In particular, these mesons may not be simple quark-ant iquark 
bound states but e.g. might contain sizable admixture of 2quark-2antiquark pairs |T^. In 
that respect we might interpret our results as a quark-antiquark model prediction for scalar 
mesons. Our results suggest that such a picture is too simple for these mesons. One might 
also speculate that the adopted ladder approximation could be insufficient. 

For the axialvector mesons we have two channels that are distinguished by their charge 
conjugation properties, c/. tables and ^ The quark-antiquark pairs that are bound to 
axialvector modes with negative charge conjugation eigenvalue tend to be lighter than those 
with the positive eigenvalue but otherwise equal quantum numbers. Generally we find that 
our predictions are lower than the assignments made by the particle data group 



We recognize from our results that the model predictions change only slightly within the 
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D 


niu 


nis 


M - 


Mu-s 


Ms-s 


0.40 


45.0 


5 X 10"^ 


0.120 


0.700 


0.917 


1.096 


0.45 


25.0 


5 X 10-3 


0.120 


0.675 


0.908 


1.099 


0.50 


16.0 


5 X 10-3 


0.115 


0.645 


0.903 


1.113 



TABLE IV. Results for the scalar mesons. The subscripts of M denote the flavor content. All 
units are GeV. 



u 


D 


■niu 


nis 


Muu 


Mus 




0.40 


45.0 


5 X 10-3 


0.120 


0.804 


0.994 


1.128 


0.45 


25.0 


5 X 10-3 


0.120 


0.858 


1.047 


1.182 


0.50 


16.0 


5 X 10-3 


0.115 


0.912 


1.098 


1.230 


experiment 1^ 


1.230 


1.270 


1.170 ? 



TABLE V. Results for the axial- vector {J^'^ = 1+-) mesons. The question mark indicates 
that the PDG did not assign the charge conjugation property of the respective resonance. All 
units are GeV. 



UJ 


D 


m„ 


nis 


Muu 


Mus 


Mss 


0.40 


45.0 


5 X 10-3 


0.120 


0.917 


1.117 


1.253 


0.45 


25.0 


5 X 10-3 


0.120 


0.918 


1.124 


1.270 


0.50 


16.0 


5 X 10-3 


0.115 


0.927 


1.140 


1.292 


experiment [Isl 


1.230 


1.270 


1.282 



TABLE VL Results for the axial-vector {J^^ = 1++) mesons. All units are GeV. 
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UJ 


D 


ruu 




Muu 


Mus 


Mss 


0.40 


45.0 


5 X 10-3 


0.120 


0.807 


0.990 


1.131 


0.45 


25.0 


5 X 10-3 


0.120 


0.861 


1.040 


1.185 


0.50 


16.0 


5 X 10-3 


0.115 


0.915 


1.085 


1.233 



TABLE VII. Results for the axial-vector mesons allowing for mixing of the Dirac structures in 
eqs. (Inl) and (pi). All units are GeV. 



large range of considered model parameters. This confirms that meson static properties are 
not too sensitive to the conjectural parameter dependence of the timelike quark propagator 
functions. Presumably meson properties whose computation involves larger timelike will 
exhibit a stronger sensitivity. 

Already from table we observe that by increasing uj the predicted mass of the J^*^ = 
1+- meson with pion flavor quantum numbers approaches the empirical mass^. We therefore 
further increased u according to the rules discussed above. For u ~ O.SGeV we reproduced 
the empirical value for the mass in that channel. However, this happened at the expense of 
significantly lowering fx and loosing the proper description of the vector mesons. We recall 
that the parameter u has a physical interpretation as the location of the maximum of the 
interaction. Thus u = O.SGeV seems intuitively too large for low-energy hadron physics 
and an unsatisfactorily description of the 0~ and 1~ mesons comes without surprise. 

For non-diagonal flavor structures such as us, charge conjugation actually is not a sen- 
sible quantum number and the corresponding axial vector mesons 1"^^ and l"*"- may mix. 
In table |V11| we present the results obtained from the full calculation that combines the 
Dirac decompositions (|ll]) and (p!2D. Since our Bethe-Salpeter formalism only yields the 
lowest mass eigenstate within a given channel, the results presented in table |V11| should be 
compared to those in table 0. The tiny changes for the flavor diagonal mesons are numerical 
artifacts. Surprisingly the changes for the non-diagonal flavor structure are also only of the 
order 1%. This suggests only a small mixing between the 1"'""'" and 1+" states with the flavor 
structure us and the l"*"^ and l"*"- channels represent good approximations to the actual 
eigenstates. 

We can extend our model beyond the light flavors up, down and strange. The only 
modiflcation is the increase of the current quark mass, mg. In Fig. ^ we show how the 
meson masses increase as the (equal) quark masses are increased into the charm sector 
rric = 1.125GeV. Exactly the same numerical code is used to construct these solutions to 
the Schwinger-Dyson and Bethe-Salpeter equations as for the light flavors. Clearly seen is 
the smooth way the masses increase from the chiral limit (mo = 0) into the heavy quark 
sector (mo = m^). This represents a convincing indicator for the stability of our technique. 
The cc-meson masses can be loosely extracted (table [VIII|) and the data are surprisingly 
well reproduced. We did not expect to be able to describe such a heavy system with such 
a simple model, derived (and fltted) as it is from pion physics. The lack of the correct UV 



^In all other channels the dependence on the parameters is much softer within the allowed window. 
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f ' ' ' ' ' 1 

0.4 0.8 1.2 

FIG. 3. Meson masses as a function of the (equal) quark mass, uiq. to = O.bGeV, D = IGGeV 
All units are GeV. 



jPiC) 




experiment |18| 


0- 


2.97 


r/c: 2.98 


1" 


3.13 


J/tp: 3.10 


0+ 


3.32 


Xco: 3.42 


1++ 


3.38 


Xci- 3.51 


1+- 


3.31 


? 



TABLE VIII. Results for the cc-meson states, lv = 0.4, D = 45.0, ttt-c = 1.125. rric is fitted 
approximately from the r/c mass. 

behavior for the gluon is seemingly at odds with the scales present. However, the present 
results suggest that the Bethe-Salpeter equation is capable of describing all the angular 
momentum states equally well in the charm quark sector. 

V. SUMMARY AND OUTLOOK 

In this paper we have studied the low-lying mesons as quark-ant iquark bound states 
in a covariant approach using an effective interaction. This interaction is characterized 
by gluon exchange with the gluon propagator being dressed by a Gaussian shape func- 
tion. The interaction is completed by the quark-gluon vertex that we take to be the tree- 
level perturbative one. In this manner the rainbow-ladder approximation to the system 
of Schwinger-Dyson and Bethe-Salpeter equation accounts for chiral symmetry. With this 
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effective interaction, we have then consistently treated this system of integral equations 
by precisely implementing the quark propagator functions that solve the Schwinger-Dyson 
equations into the Bethe-Salpeter equations. Once the effective interaction exceeds a cer- 
tain strength, the Schwinger-Dyson equations exhibit dynamical chiral symmetry breaking 
and the pseudoscalar mesons emerge as would-be Goldstone bosons. We have then used 
observed properties of the pseudoscalar mesons to determine the model parameters. The 
kaon decay constant represents a model prediction. It turned out to be in good agreement 
with the empirical data. Furthermore our results for the vector meson masses match the 
experimental data. The situation in the scalar channel is less satisfying. As we solely con- 
sider the mesons as bound states of quark-antiquark pairs, it is not surprising that the mass 
eigenvalues increase with the strangeness content. On the other hand it is astonishing that 
for current quark masses, tuq > 0.2GeV, the lightest scalar mesons turn out to be heavier 
than the lightest vector mesons. When discussing these results it must be noted that the 
role of the scalar mesons is still under intense debate. In particular, the question whether 
they should indeed be considered as quark-antiquark bound states is not yet completely re- 
solved. There are indications, see e.g. Ref. 0] and references therein, that the scalar meson 
masses should actually decrease with the strangeness content of these mesons. This can be 
understood if these mesons are considered as 2quark-2antiquark bound states in the sense 
of diquark-antidiquark systems [1^]. 

As an outlook we mention that there is an elegant way to extend the present model to 
incorporate such degrees of freedom. The Bethe-Salpeter treatment can be straightforwardly 
extended to study bound states of diquark-antidiquark pairs, once a binding mechanism is 
established. This could either be achieved by a gluon exchange similar to eq. (|1]) or by 
quark exchange between a quark and a diquark. The latter approach has been intensively 
studied and the corresponding vertex is known from modeling baryon properties [|TU|. It 
will also be interesting to see whether these additional degrees of freedom will also affect 
the mass predictions for the axialvector mesons that currently tend to be on the low side. 
Investigations in this direction are in progress. 

Finally we would like to repeat that the stability of our treatment - measured by the 
actual momentum partition invariance - allows us to even describe mesons with charm 
quantum numbers. Without further modifications of the model parameter, except the cor- 
responding current quark mass, our approach reproduces the mass eigenvalues of the flavor 
neutral states unexpectedly well. 
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FIG. 4. Dependence of the pion (left) and kaon (right) masses on the momentum partition ^. 
Model parameters are: oj = 0.4, D = 45, = 5 x 10^"^ and = 0.120. All units are GeV. For 
the calculation of the pion mass A'^ch = 4 has been used. 



APPENDIX A: RELATIVISTIC COVARIANCE 



As mentioned in section III the truncation of taking A^^ch finite in eq. ([T3| ) violates rela- 
tivistic covariance. Studying the dependence of our results on the momentum partitioning 
parameter ^ we can provide a measure of these violations originating from the numerical 
method. The results for the pseudoscalar case are shown in Fig. ^. The pion mass is ex- 
tremely stable over the whole range ^ G [0, 1], whereas for the kaon, the asymmetry between 
the quark masses must be compensated for by increasing the number of Chebyshev moments, 
Ach- Clearly, for the light pseudoscalars, the technique of inserting the quark propagators 
directly results in extremely stable results. In the case of the other mesons, the stability 
becomes less evident as one is studying heavier mass resonances (but with the same quarks). 
As an example, we display the axial-vector {J^^ = us meson (Fig. There is still 
a clear range of stability ^ G [0.4,0.7]. Outside this range however, we notice that the 
technique does not simply break down (a mass can still be found) - it merely gets numer- 
ically more difficult to get trustworthy results. We assert the increasing dependence on ^ 
with larger meson masses to the fact that with larger meson masses the calculation becomes 
more sensitive to the violent behavior seen in the timelike quark propagator functions. The 
larger the influence of that behavior to the calculation the higher we have to assume the 
number of Chebyshev moments for an adequate representation. Away from this timelike 
regime the quark propagator functions are smooth. This results in a good representation of 
the Bethe-Salpeter kernel and vertex functions with only a few Chebyshev moments. 
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FIG. 5. Dependence of the axial-vector (J^*^ = 1^ ) mass on the momentum partition One 
set of parameters has been used: lo = 0.4, D = 45, m„ = 5 x 10^'^ and = 0.120. Ah units 
are GeV. 
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